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Abstract 

Charged particle momentum distributions are studied in the reaction e + e~ — * hadrons, using 
data collected with the OPAL detector at centre-of-mass energies from 192 GeV to 209 GeV. 
The data correspond to an average centre-of-mass energy of 201.7 GeV and a total integrated 
luminosity of 433 pb _1 . The measured distributions and derived quantities, in combination 
with corresponding results obtained at lower centre-of-mass energies, are compared to QCD 
predictions in various theoretical approaches to study the energy dependence of the strong 
interaction and to test QCD as the theory describing it. In general, a good agreement is found 
between the measurements and the corresponding QCD predictions. 
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1 Introduction 



Charged particle momentum distributions are studied in the hadronic decays of virtual pho- 
tons or Z bosons (hereafter referred to as (Z/7)* decays) at centre-of-mass (cm.) energies, 
between 192 and 209 GeV, the highest available e + e~ energies. These distributions can be mea- 
sured with relatively high precision. Combined with lower energy data they provide a powerful 
test of quantum chromo dynamics (QCD) as the theory describing the strong interaction. 

Previous studies using e + e _ annihilation data at cm. energies from 130 up to 189 GeV have 
shown that QCD based models and calculations give a good description of most of the measured 
distributions @,|,|,|,||,|]. With the yet higher energy data used in this paper, we obtain even 
more stringent tests of the models and theory. In addition we use these highest energy data in 
combination with earlier measurements at lower cm. energies to study the energy dependence 
of the strong interaction. 

We measure the charged particle momentum distribution, 1/cr • d<7 c h/dx p , and the £ p dis- 
tribution, 1/cr • dc c h/d£p, where x p = 2p/y / i, £ p = ln(l/x p ), p is the particle momentum, a is 
the cross-section for non-radiative hadronic (Z/7)* decays and <r c h is the charged particle cross- 
section in these events. We also measure the distribution of the rapidity, y = ^ In ( ^pjj ) > 
where pn is the momentum component parallel to the thrust axis and E is the energy of the 
particle, and the distributions of the 3-momentum components in, p 1 ^, and perpendicular to, 
Pj^, the event plane. This plane is defined by the eigenvectors associated with the two largest 
eigenvalues of the momentum tensor, as in reference [ffl]. 

We present comparisons of the fractional momentum distributions 1/cr • dcr c h/dx p and 1/cr • 
dc c h/d^ p in the range ^/i = 14—202 GeV to QCD predictions in different theoretical approaches. 
As in earlier publications (l],[2|,|3| we also study the evolution with the cm. energy of the peak 
position £0 in the £ p distribution. 

In Section ^ a brief description of the OPAL detector is given. The samples of data and 
simulated events used in the analysis are described in Section ||[ In Section ||] the event selection 
and analysis procedure are described. In Section || we present the observed distributions and, 
together with those obtained at lower energies, compare them to the predictions of different 
Monte Carlo programs and analytic QCD predictions. A summary and conclusions are given in 
Section || 



2 The OPAL detector 



The OPAL detector was operated at the LEP e + e~ collider at CERN. A detailed description 
can be found in J?]] . The analysis presented here relies mainly on the measurement of momenta 
and directions of charged particles in the tracking chambers and of energy deposited in the 
electromagnetic calorimeters of the detector. 
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All tracking systems were located inside a solenoidal magnet which provided a uniform axial 
magnetic field of 0.435 T along the beam axis^- The magnet was surrounded by a lead glass 
electromagnetic calorimeter and a hadron calorimeter of the sampling type. Outside the hadron 
calorimeter, the detector was surrounded by a system of muon chambers. There were similar 
layers of detectors in the forward and backward endcaps. 

The tracking system consisted of a silicon microvertex detector, an inner vertex chamber, a 
large volume jet chamber, and specialised chambers at the outer radius of the jet chamber to 
improve the measurements in the z direction. The main tracking detector was the central jet 
chamber. This device was approximately 4 m long and had an outer radius of about 1.85 m. It 
had 24 sectors with radial planes of 159 sense wires spaced by 1 cm. The momenta p of tracks 
in the r-<j) plane were measured with a precision a p /p = \/0.02 2 + (0.0015 -p[GeV/c]) 2 . 

The electromagnetic calorimeters in the barrel and the endcap sections of the detector con- 
sisted of 11704 lead glass blocks with a depth of 24.6 radiation lengths in the barrel and typically 
22 radiation lengths in the endcaps. The barrel section covered the angular region | cos 6\ < 0.82 
and the endcap section the region 0.82 < | cos#| < 0.98. 



3 Data and Monte Carlo samples 



The data used in this measurement were recorded in 1999 and 2000, the final two years of the 
LEP 2 program. During this time e + e~ interactions were collected at cm. energies in the range 
192-209 GeV. The total integrated luminosity of the data sample, as evaluated using small angle 
Bhabha collisions, is 433 pb _1 . The luminosity collected at the different cm. energies is detailed 
in Figure [l] and in Table |l[ 

Monte Carlo event samples were generated at cm. energies of 192, 196, 200, 202, 206 and 
208 GeV and were processed using a full simulation of the OPAL detector |8|. Distributions 
at intermediate energies were obtained by linear interpolation between the distributions at the 
nearest neighbouring available cm. energies. In general, the Monte Carlo samples contain at 
least 10 times the statistics of the data. Signal e + e~ — ► (Z/7)* — ► qq(g) events were generated 
using the PYTHIA 6.125 f||l0| Monte Carlo program for the parton shower and hadronisation 
stage, which was interfaced with the KK2F program [|ll[] to obtain a more accurate description 
of initial state radiation (ISR). The simulation parameters were tuned to OPAL data taken at 
the Z peak p2[ . The PYTHIA Monte Carlo model has previously been shown to provide a good 
description of e + e~ annihilation data for cm. energies up to 189 GeV (see e.g. J|,|2|,||,||,[5|,[||l2|]). 

As an alternative to the string fragmentation model implemented in PYTHIA, events were 



generated using the HERWIG 5.9 [jl3[ Monte Carlo program, also tuned to OPAL data as 
described in [^,|l2|. The HERWIG Monte Carlo program implements the cluster fragmentation 
model. 

In addition, we generated events of the type e + e~ — ► 4 fermions (diagrams without interme- 
diate gluons). These 4-fermion events, in particular those with four quarks in the final state, 
constitute the major background in this analysis. Simulated 4-fermion events with hadronic and 
leptonic final states were generated using the GRC4F 2.1 Monte Carlo program. The final 

1 In the OPAL coordinate system the z axis points in the direction of the electron beam, 8 is the polar angle 
with respect to this axis and <j> is the azimuthal angle. 
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states were produced via s-channel or t-channel diagrams, including W + W and ZZ produc- 



tion. This generator was interfaced to JETSET 7.4 || using the same parameter set [12] for 
the fragmentation and decays as used for (Z/7)* events. The JETSET program implements the 
string fragmentation model. 

Two other possible sources of background events were simulated. Hadronic two-photon 



processes were evaluated using the VERMASEREN |15|, HERWIG and PHOJET |16| programs. 



Production of e + e — ► (Z/7)* — ► r + r was evaluated using the KK2F event generator in 



combination with the TAUOLA [17| decay library to simulate the decay of the tau leptons. 



In addition to PYTHIA and HERWIG we also used the ARIADNE 4.08 event gener- 
ator to compare to our corrected data distributions. The ARIADNE program implements the 
colour dipole model to describe the parton shower process. For the fragmentation stage the 
generator was interfaced to the JETSET 7.4 program. The parameter set used for ARIADNE 
is documented in [fill]. This model provides a good description of global e + e~ event properties 
at y/s = M z , as do PYTHIA and HERWIG. 



4 Data analysis 
4.1 Selection of events 

The majority of hadronic events produced at cm. energies above the Z resonance are radiative 
events in which initial state photon radiation reduces the invariant mass of the hadronic system 
to about Mz- An experimental separation between radiative and non-radiative events is therefore 
required. In addition, at the present energies above 190 GeV, ZZ and W + W~ production are 
kinematically possible and form a significant background to the (Z/7)* — ► qq events. In this 
paper we use similar techniques to those of our previous analyses of c + c~ annihilation data at 
lower cm. energies to select non-radiative (Z/7)* — ► qq events. 

4.1.1 Preselection 



Hadronic events are identified using criteria as described in [2C]. The efficiency of selecting 
non-radiative hadronic events is greater than 98% at all cm. energies, as can be seen in Table @. 
We define as particles tracks recorded in the tracking chambers and clusters recorded in the 
electromagnetic calorimeter. The tracks are required to have transverse momentum relative to 
the beam axis, px > 150 MeV/c, a number of hits in the jet chamber, iVhits > 40, a distance 
of the point of closest approach to the collision point in the r-cf) plane, do < 2 cm, and along 
the z axis, zq < 25 cm. The clusters in the electromagnetic calorimeter are required to have a 
minimum energy of 100 MeV in the barrel and 250 MeV in the endcap sections. 



4.1.2 ISR-fit selection 

To reject radiative events, we determine the effective cm. energy v 7 ? of the observed hadronic 



system as follows [21]. 
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First isolated photons are identified by looking for energy deposits greater than 3 GeV in 
the electromagnetic calorimeter, having less than 1 GeV of additional energy deposited in a 
cone of 0.2 radians around its direction. The remaining particles are formed into jets using 
the Durham p3| algorithm with a value for the resolution parameter y cu t = 0.02. A matching 
algorithm [23] is employed to reduce double counting of energy in cases where charged tracks 



point towards electromagnetic clusters. The energy of additional photons emitted close to the 
beam direction is estimated by performing three separate kinematic fits assuming zero, one or 
two such photons. Of the acceptable fits, the one with the lowest number of photons is selected. 
The value of Vs* is computed from the fitted momenta of the jets, excluding photons identified 
in the detector or close to the beam directions. The value of yfs' is set to \fs if the fit assuming 
zero initial state photons was selected. The 4-momenta of all measured particles are boosted 
into the rest frame of the observed hadronic system. 



To reject events with large initial-state radiation (ISR), we require vs' > \fs — 10 GeV. This 
is referred to as the "ISR-fit" selection. In Figure ^a we compare the distribution of the 
data set at y^^OO GeV, after the preselection was applied, to simulated (Z/7)* events and 
to 4-fermion and other background events. Simulated (Z/7)* events are classified into radiative 
events, \/ s( ruc < y/s — 1 GeV, where \J s' true is the true effective cm. energy, and non-radiative 
events, which is the complement. Though about 27% of the selected (Z/7)* events are by this 
definition radiative events, the fraction of selected events with \J s[ ruc < ,/s - 10 GeV is only 
5%. The background from 4-fermion events and the efficiency of selecting non-radiative events 
are given in Table [I]. 



4.1.3 Final selection 



The estimated background from e + e~ — > t + t~ and two-photon events of the type 77 — > qq is 
small at high yfp. To reject this background further and to ensure events are well contained in 
the OPAL detector we require at least seven accepted tracks, and the cosine of the polar angle 
of the thrust axis, | cos 9t\ < 0.9. The background from the t + t~ and two-photon events in the 
final selected sample is estimated from Monte Carlo samples to be about 0.1% and is neglected. 

To reduce the background of 4-fermion events in the remaining sample, we test the compati- 
bility of the events with QCD-like production processes. A QCD event weight Wqcd is computed 
as follows. We force each event into a four-jet configuration in the Durham jet scheme and use 



the EVENT2 [24] program to calculate the 0{o? s ) matrix element \M(pi,P2,P3,P4:)\ 2 for the 



processes e + e — » qqqq, qqgg [25], where pi, P2, P3 and P4 are the momenta of the reconstructed 



jets. Since neither quark nor gluon identification is performed on the jets, we calculate the ma- 
trix element for each permutation of the jet momenta and use the permutation with the largest 
value for the matrix element to define the event weight, 

Wqcd = , max log ( \M(pi,p 2 ,P3,P4)\ 2 ) ■ (1) 

{Pl,P2,P3,P4j V ' 

Note that the definition of the event weight contains kinematic information only and is inde- 
pendent of the value of a s . 

The weight Wqcd is expected to have large values for processes described by the QCD matrix 
element, originating from (Z/7)* — > qq, and smaller values for W + W~ or ZZ events. In Figure ||b 
we compare the data distribution of Wqcd, after the ISR-fit selection, to the expectations of our 
simulation. A good separation between the (Z/7)* and 4-fermion events is achieved by requiring 
Wqcd > -0.5. 
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In addition to the cut on Wqcd> used also in our earlier papers, we apply a likelihood based 
rejection of hadronic and semi-leptonic ZZ and W + W~ decays (see [26 and references therein). 
The distributions for the semi-leptonic and hadronic likelihoods, after all selection cuts described 
above, are shown in Figure ||a and [|b, for the data and for the signal and background Monte 
Carlo expectations. Events are rejected if the semi-leptonic likelihood is greater than 0.5 or 
the hadronic likelihood is greater than 0.25, thus reducing the remaining 4-fermion background 
by more than a factor two while reducing the expected signal by only about 2.2%. After the 
final selection criteria the estimated efficiency for selecting signal events is around 76% while 
the remaining 4-fermion background in estimated to be around 5%. The individual numbers for 
each of the six energy ranges are given Table @. 



4.2 Correction procedure 

The remaining 4-fermion background in each bin of each observable was estimated by Monte 
Carlo simulation and is subtracted from the observed bin content. A bin-by-bin multiplication 
procedure is then used to correct the observed distributions for the effects of detector resolution 
and acceptance as well as for the presence of remaining radiative (Z/7)* events. A bin- by-bin 
correction procedure is suitable for the measured distributions as the effects of finite resolution 
and acceptance cause only limited migration (and therefore correlation) between bins and the 
data are well described by the Monte Carlo programs used to determine the corrections. 

For the multiplicative correction each bin, after background subtraction, is corrected from 
the "detector level" to the "hadron level" using two samples of Monte Carlo (Z/7)* events at 
each cm. energy. The hadron level sample does not include initial state radiation or detector 
effects and allows all particles with lifetimes shorter than 3 x 10~ 10 s to decay. The detector level 
sample includes full simulation of the OPAL detector and initial state radiation and contains 
only those events which pass the same cuts as are applied to the data. The bin-by-bin correction 
factors are derived from the ratio of the distributions at the hadron level to those at the detector 
level. The correction factors tend to be around 1.15 in most bins, with values further away from 
1 at the edges of the distributions. 



4.3 Systematic uncertainties 



The experimental systematic uncertainties are estimated by repeating the analysis with varied 
experimental conditions. To reduce statistical fluctuations in the magnitudes of the systematic 
uncertainties the average over three neighbouring bins is taken. In addition, all results and errors 
obtained for the different energy ranges are combined in a single set of results, as described in 



Section 4.4, thus further reducing any statistical fluctuations in the systematic uncertainties. 



Possible inadequacies in the simulation of the response of the detector in the endcap regions 
are accounted for by restricting the analysis to the barrel region of the detector, requiring the 
thrust axis of accepted events to lie within the range |cos#t| < 0.7. This reduces the event 
sample by approximately 26%. The corresponding systematic error is the deviation of the 
results from those of the standard analysis. 

To evaluate uncertainties due to inadequacies in the track modelling the selection criteria 
for tracks were modified. The maximum allowed distance of the point of closest approach of 
a track to the collision point in the r-<j> plane, do, was changed from 2 to 5 cm, the maximum 
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distance in the z direction, zq, from 25 to 10 cm and the minimum number of hits from 40 to 
80. These modifications change the number of selected tracks by up to approximately 12%. In 
addition, track momenta in the Monte Carlo samples were smeared to degrade the resolution 
by 10%. The quadratic sum over the deviations from the standard result, obtained from each 
of these variations, is included in the total systematic uncertainty. 

Uncertainties arising from the selection of non-radiative events are estimated by repeating 
the analysis using a different technique to determine the value for y/s'. This technique differs 
from our standard algorithm in that in this case the kinematic fit always assumes one 
photon, either unobserved close to the beam direction, or detected in the ECAL. The final event 
sample with this \J~s' algorithm has an overlap of approximately 97% with the standard sample 
and is approximately the same size. The difference relative to the standard result is included in 
the total systematic error. 

Systematic uncertainties associated with the subtraction of the 4-fermion background events 
are estimated by varying the cut on Wqcd and on the 4-fermion likelihoods. We use Wqcd > 
—0.8 which increases the event sample by approximately 2%, and Wqcd > which reduces the 
event sample by approximately 8%. We also vary the cuts on the likelihood values between 0.1 
and 0.4 for hadronic 4-fermion events and between 0.25 and 0.75 for semi-leptonic events. The 
maximum deviation from the standard result obtained from the variations of the cuts on Wqcd 
and the hadronic 4-fermion likelihood is included in the total systematic uncertainty, as is the 
maximum deviation obtained from the variations of the cut on the semi-leptonic likelihood. 

In addition, we vary the predicted background to be subtracted, up and down by 5%, slightly 
more than its measured uncertainty at \fs= 189 GeV of 4% |26j| , and include the largest difference 
from the standard result in the systematic uncertainty. 

The difference in the results when we determine the bin-to-bin correction factors from sim- 
ulated (Z/7)* events generated using HERWIG instead of PYTHIA is taken as the uncertainty 
in the modelling of the (Z/7)* events. 

All systematic uncertainties are added in quadrature to obtain the total systematic errors. 
The largest contribution to the total uncertainty is that obtained using HERWIG as an alter- 
native hadronisation model in the correction procedure. Other significant contributions to the 
systematic uncertainties arise from varying the quality criteria on tracks and from varying the 
4-fermion rejection cuts. 

4.4 Combining results 

The event selection, the correction procedure and the estimate of systematic uncertainties are 
performed independently for each of the six energy regions detailed in Table [p. Afterwards the 
results are combined to obtain a single set of results with the best possible statistical precision. 
The combination procedure assumes all systematic uncertainties to be fully correlated between 
the energy regions. This means data points obtained at different energies are combined into a 
weighted average, using the statistical errors to determine the weights. 

To obtain the systematic uncertainty on the combined results all individual contributions 
to the systematic uncertainty, described in Section are first determined separately in each 
energy region. For the combined results the individual contributions to the systematic uncer- 
tainty are then taken as the average over the corresponding values obtained in each of the six 
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energy regions. Finally, the combination of the thus averaged individual contributions into the 
total systematic uncertainty on the combined results is done, as in Section <L3, by adding all 
contributions in quadrature. 



5 Results and comparison to theory 

The p±, Pj_ ut , y, x p and £ p distributions, averaged over the full cm. energy range -y/s = 192 — 
209 GeV, are tabulated in Tables [2] to || and shown in Figures ||| and ||. These distributions 
correspond to an average cm. energy of (\fs) = 201.7 GeV with a standard deviation of 
4.8 GeV. They are compared to the corresponding predictions of the PYTHIA, HERWIG and 
ARIADNE Monte Carlo programs. In general, a reasonable agreement is observed between the 
predictions and the data for the different distributions. The HERWIG Monte Carlo program 
agrees better with the data than the PYTHIA and ARIADNE models. The p^ ut distribution is 
harder in the data than predicted by any of the Monte Carlo programs. 

There are different theoretical approaches in which fractional energy^ spectra of final state 
hadrons can be calculated. At high x, conventional perturbation theory, utilising expansions 
in powers of a s , is applicable. Non-perturbative effects occurring at the final stage of the 
hadronisation process are not calculable but they can be factorised out and left to be determined 
experimentally. At low x, perturbation theory can also be applied, provided that logarithmically 
enhanced terms (e.g. ~ ln(l/a;)) are resummed to all orders. Remaining singular effects are 
not accounted for and as a consequence the predictions obtained are formally only valid for 
asymptotically high energies. Corrections to these asymptotic predictions can be calculated as 
an expansion in powers of yTiJ. An important precondition to the validity of these asymptotic 
predictions is that the non-perturbative hadron formation process is local, i.e. the properties of 
the hadronic final state closely follow the properties of the partons before hadronisation. This 
presumed equivalence is referred to as Local Parton-Hadron Duality (LPHD) |27||. For a recent 



review of particle production in the low x region we refer to [28|. 



In the following, we compare x p and £ p distributions over a large range of cm. energies 
to theoretical predictions applicable to either the low x or high x regions, in order to test the 
validity of these QCD approaches and the LPHD hypothesis. Since the ranges in x p or £ p 
over which the different predictions are valid depend not only on the terms included in these 
predictions, but also on those that are missing, we do not quantitatively know these ranges a 
priori. Our approach, therefore, is to fit predictions to the data in a range where they can be 
observed to describe the data well. We subsequently vary these ranges and extrapolate the fitted 
predictions outside them to establish the limits of the validity of the predictions. 



5.1 Comparisons to theory at low x 
5.1.1 Fong and Webber prediction 

At low x, QCD predicts destructive interference effects in soft gluon emissions, causing the 



asymptotic shape of the ^ distribution to be Gaussian |2Sj. The position of the peak in the 



2 The x and £ variables used in the description of the theory predictions are defined equivalently to x p and £ p , 
but refer to the fractional energies rather than the fractional momenta of particles. 
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distribution is predicted to increase linearly with r = ln(y / s/2A), where A is the QCD scale, 
defined using the one loop expression: a s (r) = 27r//?r, with (3 = 11 — 2Nf/N c , N c the number 
of colours and Nf the number of active flavours. Corrections to this asymptotic prediction were 
calculated up to 0(a s ) by Fong and Webber |pCR , yielding a skewed Gaussian shape for the £ 
spectrum: 

„ \ N(t) (k s5 (2 + k)5 2 s5 3 M 4 \ 



where S = (£ — and, for quark initiated jets: 

* = i^^/Di 1 -^^- (3) 

* - jim*(i-^(i + £\, ( 4) 

(5) 
(6) 

with p = ll + 2N f /N 3 and a; = 1 + N f /N 3 . The term £ 

is a non-perturbative offset to the £ 

distribution. 

We have fitted Equation (||) simultaneously to the £ p distribution measured in the present 
study and in previous studies by OPAL [31], and lower energy experiments [32,33|], covering 




the range y/s = 14 — 202 GeV. As the energy scale involved in coherent parton emissions, as 
described by the formalism, is expected to be low, we apply the formalism assuming three active 
flavours. The free parameters in the prediction are A and £o- The overall normalisation iV(r) 
is also not predicted in the formalism of Equation (S). We therefore vary the normalisation 
separately at each cm. energy in the fit. 

When comparing the measured £ p distributions to the predicted £ distributions, we neglect 
the effects of the finite masses of the hadrons produced in the final state. For this reason the 
highest £ p values, corresponding to momenta lower than 300 MeV/c, are excluded from the 
fit, avoiding the region where such mass effects could play a significant role. Moreover, as the 
formalism of Equation (||) is not expected to be valid at low £ (high x), where logarithms of the 
form ln(l — x) are important, data below £ p = 2.0 are also not considered in the fit. 

Figure ||] shows the £ p data compared to the fitted predictions of Equation (||) . For clarity, 
not all OPAL data included in the fit are shown in the figure, specifically our data at 161, 172, 
183 and 189 GeV. The data in the range included in the fit (full lines) are well described by the 
prediction. At low £ p the prediction clearly lies above the data, while at high £ p it appears to 
give a reasonable description of the measurements beyond the fitted range. The total \ 2 of the 
fit is 121 for 174 degrees of freedom when including the full experimental errors (statistical and 
systematic errors added in quadrature) in the fit and assuming them to be uncorrelated. This is 
the only way to treat all data consistently since for the lower energy results only the combined 
statistical and systematic uncertainties are available. As a consequence, the quoted x 2 value 
presumably underestimates the true \ 2 °f tli e fit. 

The obtained values for A and £o are 153±7±55 MeV and —0.70 ± 0.03 ± 0.23, respectively, 
where the first errors are the fit uncertainties, defined as the square root of the diagonal elements 
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of the covariance matrix, and the second errors were obtained by varying the fit range. The 
lower limit on £ p was varied between 1.5 and 2.5 and the upper limit between the £ p values 
corresponding to particle momenta of 200 and 400 MeV/c. The errors on the parameters do not 
include an estimate of how they may be affected by missing terms in the predictions. Therefore 
these parameters are only meaningful within the formalism. A more detailed test of the formal- 
ism of Equation @ than the one presented here would need to account for more subtle effects 
such as e.g. the changing flavour composition with cm. energy. For this, additional information 
would need to be included in the fit, which goes beyond the scope of our analysis. 



5.1.2 MLLA prediction 



A prediction for the £ distribution can also be formulated in the so-called modified leading-log- 
approximation (MLLA, for a review see e.g. ]pq1). This approximation constitutes a complete 
resummation of single and double logarithmic terms Q . In the MLLA approach a so-called 
limiting spectrum prediction for the £ distribution can be formulated in which the QCD scale, 
A, and the scale at which the perturbative process terminates, Qo, are taken to have the same 
value. In this study we use the equations given in fl28| : 



^ = 2K^D^,r), (7) 
a d£ N c 

where Cp = (N% — 1) /2N C and K h is a hadronisation constant which accounts for the number of 
hadrons of type h produced per final state parton. In some earlier publications the factor of 2 in 
Equation (0), accounting for the two quark jets present in e + e~ annihilation, was incorporated 
into K h . From the LPHD hypothesis K h is expected to be independent of the underlying 
process, including its scale. For the expression for the limiting spectrum distribution, D lim , 



we refer to [28]. In this formalism r = In (y / i/2A e ff ), where A e g is an effective scale factor 
representing both the QCD scale, A, and the cut-off scale, Qq. Contrary to the Fong and 
Webber prediction, the MLLA limiting spectrum formalism predicts not only the shape of the 
£ distribution but, apart from the hadronisation correction, also its normalisation. 

In Figure [?]a, the same £ p data shown in Figure ^ are compared to the fitted predictions of 
Equation (0), where the free parameters are the effective scale, A e g , and a separate hadronisation 
constant, K h , for each cm. energy. We choose to fit a separate hadronisation constant at 
every energy in order subsequently to examine the obtained K h values to establish how well 
the formalism and the LPHD hypothesis work. Higher order effects, not accounted for in the 
formalism, are expected to give rise to some energy dependence in the K h factors, in particular 



at the lowest cm. energies (see e.g. [35]). As in Section 5.1.1 , we apply the formalism assuming 
three active flavours. The fit is restricted to the region around the peaks defined by 0.75 + 
0.331og(^) < i p < 0.9 + 0.81og(Vs). 

We find the fit describes the data well within the fitted range (full line) as illustrated by 
the total x 2 of the fit of 122 for 132 de grees of freedom. Outside the fit range, both at low 
and high £ p , deviations of the predictions from the data are observed. The MLLA fit is seen to 
describe the low £ p data better, but the high £ p data worse, than the fit of the Fong and Webber 
prediction discussed in Section |5.1.1| (see Figure The value of A e fj obtained from the fit is 
254 ± 3 ± 31 MeV, where the first error is the fit uncertainty and the second error is obtained 
from varying the fit range. Both the lower and upper limits on £ p were lowered and raised by 
0.5 units in £ p . 
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In Figure 0b the obtained K h factors are shown as a function of the cm. energy. The errors 
shown are the fit error and the combined error from the fit and from varying the fit range, as 
above. Above 130 GeV the K h factors obtained from our fit are consistent with being constant. 
Below 50 GeV a modest rise in the values of the K h factors is observed. 



5.1.3 The energy dependence of £o 



The MLLA also provides a prediction for the energy evolution of the peak position, £o> of the £ 
distribution p3,|35|: 



1 C C *yr 

- + J +0T 2 

2 V T T 



(8) 



with t defined as in Section |5.1.1| and C = a A /WN c b 



To determine the peak position £o of the measured £ p distributions in each of the energy 
regions of the present study, we fitted Equation @ to these distributions in the range 2.0 < £ p < 
6.2, taking the position of the maximum of the fitted function as the result for £q- The systematic 
uncertainties are determined by repeating the fit using the systematic variations described in 



Section |L3|. In addition, the uncertainty due to the choice of the fit range is estimated by 
repeating the fit using two alternative fit ranges: 3.2 < £ p < 4.8 and 1.6 < £ p < 6.6, taking 
the larger deviation from the nominal result as the systematic uncertainty. The uncertainties 
are added in quadrature to obtain the total systematic uncertainty. As for the momentum 
distributions, we combine the £o results from the different energies to obtain: 

fo(201.7 GeV) = 4.158 ± 0.007(stat.) ± 0.036(syst.). 



At cm. energies below 91 GeV, with the exception of the TOPAZ result [j33| ], no £o deter- 
minations were published. Therefore, we determined £o values for all energies below 91 GeV by 
fitting Equation (||) to the published £ p (or x p ) results in a narrow region around the peak. This 
procedure ensures that the £o results are obtained in a uniform manner over the entire range 
of energies studied. The error on £o for the lower energies is taken from the variation in the 
peak position when modifying the fit range. For energies above 91 GeV this was found to be 
the dominant source of uncertainty. 

In Figure [8| our £o result is shown together with earlier £o results from OPAL |3T], |l|, ||, ||] 
and other LEP experiments [|,37,38] and the £o values we determined for the lower energy 
experiments p2 , 33 , 39 , 4C , O] . The results are compared to the predictions of PYTHIA, HERWIG 



and ARIADNE and to the fitted prediction of Equation (||). For the fit three active flavours 
were assumed and the scale, A, was the only free parameter. The fit yields A = 203 ± 2 MeV, 
where the error is the fit uncertainty. The fitted MLLA prediction is in good agreement with the 
data. The PYTHIA and ARIADNE predictions are slightly below the data, while the HERWIG 
prediction, in particular at the highest cm. energies, is much lower than the data. 



5.2 Comparisons to theory at high x 

Another theoretical framework for momentum distributions is formulated in terms of the scale 
evolution of fragmentation functions. In this framework, a factorisation between the pertur- 
bative scattering process on one hand, and the parton emission cascade and non-perturbative 
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hadronisation stage on the other hand, is formulated as follows (see e.g. [j42 , 43 , 44f ) : 
l<W ' t- E C f (z,a s ,VS)D h f (*,n). 

J x Z , , , Z 



a dx 



(9) 



f=g,u,d,s,c,b 



The Cf(x,a s , ^/s) are coefficient functions, describing the probability to obtain a parton / 
(gluon, g, or quark, u, d, s, c,b) with energy fraction x in an e + e~ collision, and Z?j(x,/x) are 
fragmentation functions, describing the expected distribution of final state hadrons h from an 
initial parton / at scale [i. The x dependence of fragmentation functions is not known a priori, 
but their dependence on the energy scale is predicted by theory. In QCD this scale dependence, 



also known as scaling violation, is described by parton evolution equations [45], or "DGLAP" 
equations. These equations, as well as the coefficient functions of Equation d9T), are currently 



known to next-to- leading order (NLO) accuracy in a s [46, 42, 43, 44] 



Using a computer program [47] written for scaling violation analyses of proton structure 
functions in NLO QCD, to which we added NLO evolution equations for fragmentation functions 
and NLO coefficient functions for e + e~ annihilation, we performed a fit of Equation (||) to 
OPAL and lower energy data. These data include the inclusive charged particle x p distributions 
measured here and at lower energies [50,|,|2|,|,3^, 39|, 48|, 49], flavour tagged x p distributions [|5C|| , 



and results on the helicity components in charged particle production [51]. 



The fit involves defining fragmentation functions Dj(x,[io) at an arbitrarily chosen input 
scale, fio, evolving these over the entire relevant range of scales using the parton evolution 
equations, and using them to calculate the theory expectation for the measured distributions 
(Equation (||)). The a priori unknown input consists of the input fragmentation functions and 
the strong coupling constant a s . The fragmentation functions are parameterised, similarly to 



reference [51], as 



xD? 



Ni(l 



-Ci (In a;) 2 



(10) 



where the parameters iVj, a«, j3j and q are determined independently for i = g,u,(ds),c,b. 
The separation between the fragmentation functions of the different quark flavours relies on the 
(uds), c and b flavour tagged data included in the fit and on the fact that as a consequence 
of the different electric and weak charges the relative contributions from d, s and b quarks 
compared those from u and c quarks, vary with the cm. energy. No data included in the fit 
provide information to separate the contributions from d and s quarks. Therefore only one 
fragmentation function is defined for these two flavours. 



In earlier studies 42 , 44 , 52, 53, 54| , similar fits were performed, including data up to cm. 
energies of 91 GeV, to obtain measurements of a s . The higher energy LEP 2 data do not add 
significantly to the achievable precision on a s , which is mostly constrained by the very precise 
low energy and LEP 1 results. We therefore choose, in our default fit, to fix a s at its current 
best estimate of a s (Mz) = 0.1181 ± 0.002 foil , to test how well the fitted predictions describe 
the data over the now extended range of scales. The data included in the fit are restricted to 
x p values between 0.07 and 0.80. In addition, the particle momenta were required to be greater 
than 1 GeV/c, effectively raising the lower x p limit for the lowest energy experiments. This cut 
is motivated by the fact that the present formalism does not account for destructive interference 
effects in soft gluon emissions (see Section [O]). 

In Figure [| x p distributions measured at cm. energies from 14 to 202 GeV are compared to 
the fitted predictions. The theory predictions are in good agreement with the data in the fitted 
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range (full lines). The total x 2 of the fit is 198 for 188 degrees of freedom. As for the fits described 
in Section |5.1| , this x 2 result is based on the full (statistical and systematic) uncertainties on the 
data and neglects any correlations between data points. Outside the fitted range, at low x p and 
low cm. energies, charged particle production is suppressed. This suppression is not described 
by the predictions, which do not include soft gluon interference effects, as stated above. 



Scaling violations in charged particle production are shown in Figure 10, where the l/er • 
da c h/dx p results are shown as a function of -y/s for a number of x p bins. Negative scaling 
violations are visible for the higher x p values. At low x p the charged particle production rate 
demonstrates little dependence on the cm. energy. Changes in the flavour composition with 
cm. energy also influence the observed level of scaling violations, especially at high x p and 
around yfs = 91 GeV where a relatively large fraction of events is expected to contain b quarks. 

The data in Figure [H] are also compared to the fitted NLO predictions. Overall, a good 
agreement between data and theory is found in the region of the fit (full lines). We thus 
conclude that inclusive charged particle production data can be described by the QCD based 
parton evolution formalism, using the world-average value for a s . When we leave a s as a 
free parameter in the fit, we obtain a s (Mz) = 0.113 ± 0.005 ± 0.007, where the first error 
is the fit uncertainty and the second the uncertainty obtained by independently varying the 
renormalisation and factorisation scales up and down by a factor 2. For the dependence of 
Equation @ and the parton evolution equations on these scales we refer to [45, 4(1 42, 43, 44]. 



6 Summary and conclusion 



We have presented a measurement of charged particle momentum distributions in hadronic e + e~ 
annihilation events produced at LEP 2 at cm. energies between 192 and 209 GeV. The results 
are combined into a single set of distributions, corresponding to an average cm. energy of 
= 201.7 GeV. 

The results, in combination with those obtained at lower cm. energies, were compared to 
QCD Monte Carlo models and to analytical QCD predictions calculated in various approaches. 
In general, good agreement is observed between the data and these theory predictions in the 
regions where the latter are expected to be valid. 
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yfs range (GeV) 


191 - 193 


195 - 197 


199 - 201 


<VS> (GeV) 


191.6 


195.5 


199.5 


Integrated luminosity (pb x ) 


29.15 ±0.09 


72.57 ±0.19 


74.51 ±0.20 


Preselection 


2796 


6820 


6262 


Expected 


2769 ± 8 


6448 ± 16 


6284 ± 15 


4-fermion background (%) 


20.1 ±0.1 


21.9 ±0.1 


24.8 ±0.1 


Eff. non-rad. events (%) 


98.4 ±0.1 


98.3 ±0.1 


98.5 ±0.1 


"ISR-fit" Selection 


765 


1848 


1734 


Expected 


745 ±4 


1793 ± 8 


1771 ± 8 


4-fermion background (%) 


30.2 ±0.3 


31.8 ±0.3 


33.6 ±0.3 


Eff. non-rad. events (%) 


89.9 ±0.3 


90.0 ±0.2 


89.8 ± 0.2 


Final selection 


507 


1111 


1047 


Expected 


461 ±4 


1074 ± 6 


1039 ± 6 


4-fermion background (%) 


4.6 ±0.1 


4.6 ±0.1 


5.3 ±0.1 


Eff. non-rad. events (%) 


77.0 ± 0.4 


76.3 ± 0.3 


76.1 ±0.3 



y/s range (GeV) 


201 - 202.5 


202.5 - 205.5 


205.5 - 209.5 


(Vi) (GeV) 


201.6 


204.9 


206.6 


Integrated luminosity (pb x ) 


37.97 ±0.11 


81.96 ±0.19 


137.14 ±0.29 


Preselection 


2971 


6031 


9847 


Expected 


3021 ± 8 


6240 ± 16 


10114 ±26 


4-fermion background (%) 


24.4 ±0.1 


25.8 ±0.2 


26.5 ±0.2 


Eff. non-rad. events (%) 


98.3 ±0.1 


98.5 ±0.1 


98.4 ±0.1 


"ISR-fit" Selection 


838 


1781 


2840 


Expected 


881 ±4 


1857 ± 8 


3052 ± 14 


4-fermion background (%) 


33.7 ±0.3 


34.5 ± 0.3 


34.8 ± 0.3 


Eff. non-rad. events (%) 


89.4 ±0.2 


89.8 ±0.2 


89.7 ±0.2 


Final selection 


494 


1091 


1647 


Expected 


514 ±3 


1075 ± 6 


1758 ± 10 


4-fermion background (%) 


5.2 ±0.1 


5.6 ±0.2 


5.7 ±0.2 


Eff. non-rad. events (%) 


75.6 ± 0.3 


75.8 ±0.3 


75.9 ±0.3 



Table 1: Data samples corresponding to six ranges in cm. energy. The luminosities and the 
luminosity weighted mean cm. energies are given in the first two rows. The rows labelled as 
"Preselection", "ISR-fit selection" and "Final selection" correspond to the number of events 
passing these criteria. For each of these selection stages the expected number of events, the 
remaining fraction of 4-fermion background and the signal efficiency are also given. The errors 
on the integrated luminosity include statistical and systematic uncertainties added in quadrature. 
All other errors are the statistical errors only. 
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± 
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± 
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± 


0.006 


± 
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± 
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± 
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Table 2: Measured values of the distribution l/a ■ d<r ch /dp^ at (y/s) = 201.7 GeV. The first 
error is statistical, the second systematic. 
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± 
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± 
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± 
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± 
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± 
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± 
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± 
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± 


0.026 



Table 3: Measured values of the distribution l/a • dcr ch /dpf * at (y/s) = 201.7 GeV. The first 
error is statistical, the second systematic. 
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Table 4: Measured values of the distribution 1/a ■ da^/dy at (y/s) = 201.7 GeV. The first 
error is statistical, the second systematic. 



Xp 


1/a 
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± 
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± 
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Table 5: Measured values of the distribution 1/a • da c h/dx p at {y/s) = 201.7 GeV. The first 
error is statistical, the second systematic. 
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Table 6: Measured values of the distribution 1/a ■ da^/d^p at (yfs) = 201.7 GeV. The first 
error is statistical, the second systematic. 
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Figure 1: Integrated luminosity collected by OPAL in 1999 and 2000 at different cm. energies. 
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QCD 

Figure 2: Distributions of (a) the effective cm. energy \fs i and (b) the QCD event weight 
Wqcd in the data at ^fs = 200 GeV (full points). The data are shown with statistical errors 
only. Estimated signal and background (BG) contributions are shown as detailed in the upper 
figure. The vertical lines indicate where the selection cuts are applied. 
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Figure 3: Distributions of the likelihood for events to be (a) a semi-leptonic 4-fermion event or 
(b) a hadronic 4-fermion event, in the data at ^fs = 200 GeV (full points). The data are shown 
with statistical errors only. Estimated signal and background (BG) contributions are shown as 
detailed in the upper figure. The vertical lines indicate where the selection cuts are applied. 
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Figure 4: Distributions of the charged particle momenta (a) p± and (b) p^ ut and of the rapidity 
(c) y&t{y/s) = 201.7 GeV, compared to the corresponding PYTHIA, HERWIG and ARIADNE 
predictions. The statistical errors and the combined statistical and systematic errors are shown 
as the inner and outer error bars, respectively. 
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Figure 5: Distributions of the charged particle fractional momenta (a) x p and (b) £ p at (y/s) = 
201.7 GeV, compared to the corresponding PYTHIA, HERWIG and ARIADNE predictions. 
The statistical errors and the combined statistical and systematic errors are shown as the inner 
and outer error bars, respectively. 
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Figure 6: £ p spectra measured in the range yfs = 14 — 202 GeV, compared to the Fong and 
Webber (3(| predictions (Equation (||) in this paper), fitted to the data. The full lines indicate 
the region of the fit. For clarity, not all OPAL data included in the fit are shown in the figure. 
The error bars on the data represent the combined statistical and systematic uncertainties. 
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Figure 7: a) £ p spectra measured in the range ^/s = 14 — 202 GeV, compared to MLLA 
predictions in the limiting spectrum approach ( p8| , Equation (0) in this paper), fitted to the 
data. The full lines indicate the region of the fit. For clarity, not all OPAL data included in the 
fit are shown in the figure. The error bars on the data represent the combined statistical and 
systematic uncertainties. 

b) The hadronisation constants K h at different y/s as obtained from the fit, described in Sec- 
tion 5.1.2 . The inner error bar represents the fit error and the outer error bar the combined fit 
error and the error obtained from varying the fit range. 
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Figure 8: Evolution of the position of the peak in the £ p distribution, £q> with cm. energy, in 
the range y/s = 14 — 202 GeV, compared to the fitted MLLA prediction (Equation (||) in this 
paper) and to predictions of PYTHIA, HERWIG and ARIADNE. The error bars of the data 
represent the combined statistical and systematic uncertainties. 
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Figure 9: x p spectra measured in the range y/s = 14 — 202 GeV, compared to fitted NLO 
predictions (Equation (^) in this paper) using the world average value of a s (Mz) = 0.1181. The 
full lines indicate the region of the fit. The uncertainties of the data represent the combined 
statistical and systematic errors. 
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Figure 10: Inclusive charged particle rate, 1/a • da c h/dx p , as a function of y/s in bins of x p , 
compared to fitted NLO predictions (Equation (^) in this paper) using the world average value 
of q s (Mz) = 0.1181. The full lines indicate the region of the fit. The uncertainties of the data 
represent the combined statistical and systematic errors. 
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